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Abstract 

A  simple  two-dimensional  model  involving  kinetics  and  mass  transfer  in  a  PEM  fuel  cell  cathode  is  presented.  In  this  model,  the  catalyst 
layer  was  simplified  as  an  infinite  thin  film.  The  oxygen  mass  transfer  in  the  gas  diffusion  layer  (GDL)  was  described  using  a  pure  diffusion 
equation  that  introduced  equivalent  oxygen  diffusivity.  The  PEM  fuel  cell  performance  under  the  influence  of  current  collector  ribs  was 
investigated.  The  results  show  that,  the  existence  of  ribs  causes  the  GDL  to  be  used  only  partly  in  the  mass  transfer  process.  The  GDL 
effectiveness  decreases  with  the  cell  current  density  and  increases  with  the  width  of  the  gas  flow  channels.  The  PEM  fuel  cell  performance 
decreases  with  an  increase  in  GDL  thickness  if  the  GDL  porosity  is  low.  However,  when  a  high-porosity  GDL  is  used,  the  optimal  thickness 
becomes  an  indicator  determining  the  maximal  PEM  fuel  cell  performance. 
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1.  Introduction 

Most  PEM  fuel  cell  designs  are  based  on  planar  and  repet¬ 
itively  stacked  structures.  Each  cell  in  the  stack  has  two  bipo¬ 
lar  plates  pressed  against  the  membrane  electrode  assembly 
(ME  A),  as  shown  in  Fig.  1.  The  ME  A  is  the  core  component 
of  PEM  fuel  cells,  which  consists  of  the  proton  exchange 
membrane  (PEM),  anode  and  cathode  electrodes.  An  elec¬ 
trode  comprises  both  the  gas  diffusion  layer  (GDL)  and  the 
catalyst  layer.  The  GDL  is  usually  made  of  carbon  papers  or 
carbon  cloths,  forming  the  outmost  portion  of  the  ME  A  and 
positioned  next  to  the  bipolar  plate.  The  gas  flow  channels 
grooved  on  both  sides  of  the  bipolar  plate  are  designed  to 
distribute  the  reactant  gas  to  the  electrode  reaction  sites.  The 
current  collector  ribs  laid  between  two  neighbouring  chan¬ 
nels  are  the  paths  for  the  electric  current. 

Because  the  electrode  portion  covered  by  the  ribs  is  not 
directly  exposed  to  the  channels,  it  suffers  from  a  slow  reac¬ 
tant  gas  mass  transfer.  The  ribs  can  be  regarded  as  barriers 
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to  mass  transfer;  however,  they  are  indispensable  to  electric 
current  conduction  in  a  PEM  fuel  cell. 

A  complete  understanding  of  the  mass  transfer  phenom¬ 
ena  within  the  GDL,  under  the  influence  of  current  collector 
ribs,  will  facilitate  a  proper  PEM  fuel  cell  design.  However, 
there  have  been  very  few  studies  in  the  literature  dealing  with 
such  a  critical  problem.  West  and  Fuller  [1]  studied  the  ef¬ 
fects  of  rib  sizing  and  the  GDL  thickness  on  the  current  and 
water  distributions  within  a  PEM  fuel  cell.  They  found  that 
the  ribs  only  slightly  altered  the  cathode  potential  for  a  given 
current  density,  but  had  a  significant  influence  on  water  man¬ 
agement.  Hental  et  al.  [2]  experimentally  investigated  the 
effects  of  both  rib  and  channel  widths  on  the  performance 
of  single  PEM  fuel  cells.  Recently,  Yan  et  al.  [3]  developed 
a  two-dimensional  mass  transport  model  to  investigate  the 
anode  gas  flow  channel  cross  section  and  GDL  porosity  ef¬ 
fects.  They  found  that  an  increase  in  either  the  GDL  porosity, 
channel  width  fraction,  or  the  number  of  channels  could  lead 
to  better  cell  performance. 

The  PEM  fuel  cell  anode  overpotential  is  negligible  in 
comparison  with  the  cathode  overpotential,  thus  the  hydro¬ 
gen  mass  transport  at  the  anode  is  considered  less  significant 
than  the  oxygen  mass  transport  at  the  cathode.  The  purpose  of 


42 


K.T.  Jeng  et  al.  /  Journal  of  Power  Sources  138  (2004)  41-50 


Gas  diffusion  Catalyst 

layer  i aver 


MEA 


Fig.  1.  Schematic  diagram  of  the  PEM  fuel  cell  and  the  enlarged  gas  diffusion  layer  element. 


this  study  is  to  investigate  the  effect  of  current  collector  ribs 
on  the  oxygen  mass  transfer  through  the  GDL.  We  will  deter¬ 
mine  the  equivalent  oxygen  diffusivity  within  the  GDL  in  an 
operating  PEM  fuel  cell  using  a  one-dimensional  model  first. 
The  so  obtained  equivalent  diffusivity  will  then  be  used  to 
describe  the  oxygen  mass  transfer  in  the  cathode  gas  mixture 
under  the  influence  of  ribs  using  a  two-dimensional  model. 
Figures  on  the  distribution  of  oxygen  concentration  within 
the  GDL  and  the  variation  in  current  density  along  the  cata¬ 
lyst  layer  will  be  recorded.  The  rib  effects  on  the  cell  perfor¬ 
mance  and  on  the  GDL  effectiveness  are  to  be  investigated. 
This  study  provides  a  novel  method  in  enhancing  the  oxygen 
mass  transfer  through  the  GDL. 


2.  Equivalent  oxygen  diffusivity  in  cathode  gas 
mixture 

In  this  study,  a  two-dimensional  model  for  oxygen  mass 
transport  within  a  GDL  has  been  simplified  into  a  pure  dif¬ 
fusion  problem  by  the  introduction  of  an  equivalent  oxygen 
diffusivity.  To  obtain  the  equivalent  oxygen  diffusivity,  the 
mass  transfer  within  the  GDL  is  investigated  using  a  one¬ 
dimensional  model  as  the  starting  point. 


2.7.  One -dimensional  model 

The  one-dimensional  GDL  model  is  schematically  shown 
in  Fig.  2.  Let  z  =  0  specify  the  GDL/channel  interface  and  z 
increase  in  the  oxygen  transport  direction.  The  gas  phase  in 
the  GDL  is  a  mixture  of  oxygen,  water  vapour  and  nitrogen 
(if  air  is  used  as  the  cathode  gas).  We  assume  that  the  cathode 
gas  mixture  acts  as  an  ideal  gas  and  the  fuel  cell  operates 
under  steady- state  conditions.  The  Stefan-Max  well  equation 
for  multi-component  diffusion  can  be  used  to  describe  the 


gradient  in  the  component  mole  fractions: 
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where  Xi  and  N[  represent,  respectively,  the  mole  fraction  and 
the  molar  flow  rate  of  species  i.  Df?,  the  effective  binary 
diffusivity  of  the  gas  pair  i-j  in  the  porous  media,  and  Pu  the 
pressure  of  the  gas  mixture  in  the  GDL  which  is  considered 
to  be  constant  and  equal  to  that  in  the  gas  flow  channel.  The 
effective  binary  diffusivity  can  be  evaluated  from  the  bulk 
binary  diffusivity  D\j  using  Bruggeman’s  correction  [4]: 


Df  =  e3'2Dij 


(2) 


where  s  is  the  porosity  of  the  GDL. 


Fig.  2.  Schematic  diagram  of  the  one-dimensional  model. 
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Under  steady-state  conditions  and  from  species  conser¬ 
vation,  we  have  observed  that  the  molar  flow  rates  of  the 
oxygen,  water  vapor  and  nitrogen  are  all  constant.  Discount¬ 
ing  the  crossover  of  the  reactant  gases,  the  oxygen  molar  flow 
rate  can  be  related  to  the  current  density  I  as: 


N0  = 


I 

4 F 


(3) 


Because  nitrogen  is  inert,  we  have: 


Nn  =  0 


(4) 


The  water  vapor  transport  rate  through  the  GDL  is  the  sum 
of  the  water  production  rate  at  the  catalyst  layer,  I/2F ,  and 
the  water  transport  rate  through  the  PEM,  which  is  equal 
to  the  water  flow  rate  dragged  by  protons  (electro-osmotic 
drag)  minus  the  water  back  diffusion  rate  due  to  the  water 
concentration  gradient  in  the  PEM. 


2.2.  Net  electro -osmotic  drag  coefficient 


A  net  electro-osmotic  drag  coefficient  ft  is  defined  as  the 
number  of  net  water  molecules  transported  from  the  anode 
to  the  cathode  per  proton  [5].  This  can  be  used  to  relate  the 
net  water  flow  rate  through  the  PEM  with  the  current  density 
I.  Assuming  that,  within  the  GDL,  the  water  vapor  is  trans¬ 
ported  in  a  direction  opposite  to  z,  the  flow  rate  of  the  water 
vapor  takes  the  form  of: 

=  <5> 


The  net  water  flow  rate  through  the  PEM  has  been  inves¬ 
tigated  by  using  model  simulations  [5,6]  or  experimental 
methods  [7,8].  It  has  been  found  that  ft  is  determined  by 
the  operating  conditions  (such  as,  current  density,  tempera¬ 
ture,  pressure,  stoichiometry  and  humidity  of  the  inlet  gases), 
and  the  membrane  thickness.  Choi  et  al.  [7]  reported  that 
the  value  of  /3  decreases  with  the  current  density,  but  it  is 
nearly  constant  above  a  current  density  of  200  mA  cm-2.  For 
Nafion®  115,  the  constant  is  around  0.3  under  operating  con¬ 
ditions  of  70  °C,  1  atm  H2/O2,  and  humidified  reactant  gases. 
The  constant  goes  up  to  0.4  when  the  O2  is  not  humidified. 
Janssen  and  Overvelde  [8]  reported  that,  for  Nafion®  105,  the 
constant  is  around  0.2  under  operating  conditions  of  80  °C, 
1.5  atm  H2/air,  and  humidified  H2.  The  p  value  has  been  set 
to  0.35  when  we  calculate  the  equivalent  oxygen  diffusivity 
in  the  cathode  gas  mixture. 

Combining  Eqs.  (3)-(5)  with  Eq.  (1),  one  can  obtain  the 
following  two  differential  equations  for  v0  and  xn: 
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Because  the  cathode  gas  mixture  consists  of  oxygen,  wa¬ 
ter  vapor,  and  nitrogen,  the  following  relation  for  the  mole 
fraction  holds: 


x0  4”  T  Xn  —  1  (8) 

Both  Eqs.  (6)  and  (7)  can  be  integrated  along  with  Eq.  (8), 
from  the  initial  x0  |z=o,  xw  |z=o,  and  xn  |z=o  values,  to  yield  the 
spatial  variation  in  v0,  vw  and  xn  throughout  the  GDL.  The 
x0\z=o  and  xn\z=o  values  can  easily  be  evaluated  as  long  as 
the  xw|z=o  value  is  set. 

When  solving  Eqs.  (6)-(8),  we  assume  that  the  water  vapor 
is  saturated  at  the  GDL/catalyst  layer  interface.  Thus,  the 
chosen  xw|z=o  value  must  make  the  partial  pressure  of  water 
vapor  at  z  =  satisfy  the  following  relation: 

*wl  Z=8AP,  =  PT  (9) 

where  P^at  denotes  the  saturation  pressure  of  water  vapor  at 
the  PEM  fuel  cell  operating  temperature.  This  consequently 
results  in  a  shooting  problem  in  terms  of  the  initial  xw|z=o 
value,  which  can  be  solved  using  the  appropriate  shooting 
technique  [9]. 

Fig.  3  shows  the  v0,  vw  and  xn  profiles  in  the  one¬ 
dimensional  GDL  for  three  different  current  densities  when 
air  at  1  atm  is  supplied  to  the  cathode  channel.  The  estimated 
parameters  and  properties  used  to  describe  the  GDL  one¬ 
dimensional  mass  transfer  are  given  in  Table  1.  As  can  be 
seen  from  Fig.  3,  all  of  the  mole  fraction  graphs  are  almost 
linear  within  the  GDL  for  the  current  densities  investigated. 
Other  results  (not  shown  here)  indicate  that  the  linearity  holds 
up  to  a  GDL  thickness  greater  than  1  mm  and  higher  porosity 
dose  not  break  the  linearity  either.  The  reason  why  there  is 
a  linear  variation  in  the  mole  fractions  is  because  the  GDL 
thickness  is  very  small.  The  linearity  of  the  graphs  suggests 
an  “equivalent  diffusivity”  approach  to  resolving  the  GDL 
mass  transfer  problem. 

2.3.  Equivalent  oxygen  diffusivity  estimation 

Fick’s  law  in  terms  of  the  molar  flow  rate  of  oxygen,  as 
shown  below,  can  be  used  to  describe  the  transport  of  oxygen 
in  the  GDL: 

N0  =  X0  J2  NJ  ~  (£3/2A,-mix)^  (10) 

7=0, n,w  ^ 


Table  1 

Physical  parameter  and  property  values  used  in  one-dimensional  model 


Physical  parameter/property  Value 

Cell  temperature,  T  (K)  353.15 

Thickness  of  gas  diffusion  layer,  $d  (cm)  0.03 

Gas-pair  pressure-diffusivity  product,  PtDn-0  (atm  cm2  s-1)  0.279  [10] 

Gas-pair  pressure-diffusivity  product,  PtDw~0  (atm  cm2  s-1)  0.37  [10] 

Gas-pair  pressure-diffusivity  product,  PtDw-n  (atm  cm2  s-1)  0.387  [10] 

Net  electro-osmotic  drag  coefficient,  f3  0.35 

Porosity  of  gas  diffusion  layer,  s  0.  3 

Nitrogen-oxygen  mole  ratio  in  gas  flow  channel  3.76 
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Dimensionless  Position  (z/8d) 

Fig.  3.  Spatial  variations  of  mole  fraction  within  the  one-dimensional  GDL. 


where  ZVmix  is  the  diffusivity  of  oxygen  in  the  cathode  gas 
mixture  and  C  is  the  oxygen  concentration.  The  first  term 
on  the  right  hand  side  of  Eq.  (10)  is  the  oxygen  molar  flow 
rate  resulting  from  the  gas  mixture  bulk  motion  (convective 
transport).  The  second  term  denotes  that  resulting  from  the 
diffusion  superimposed  on  the  bulk  flow  (diffusive  transport). 
The  oxygen  diffusivity,  D0_ mix>  in  a  multicomponent  gas  mix¬ 
ture  with  a  bulk  motion  like  that  in  the  cathode  is  difficult  to 
calculate  [1 1,12],  so  we  refer  to  the  findings  in  Fig.  3  and  the 
approach  is  discussed  next. 

From  the  linear  variation  in  the  oxygen  concentration 
within  the  GDF,  we  may  simplify  the  Eq.  (10)  expression 
by  introducing  an  equivalent  oxygen  diffusivity,  D^_mix,  and 
relating  the  oxygen  molar  flow  rate  to  only  its  concentration 
gradient  as: 
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Fig.  4.  Equivalent  oxygen  diffusivity  in  the  cathode  gas  mixture. 


where,  A C  —  (x0\z=$d  —  ^0\z=o)Pt/  RT,  is  the  drop  in  oxy¬ 
gen  concentration  across  the  GDF.  The  equivalent  oxygen 
diffusivity  in  the  cathode  gas  mixture  can  then  be  calculated 
from  the  relation  between  molar  flow  rate  and  oxygen  con¬ 
centration  gradient. 

Fig.  4  shows  the  £)®_mix  results  calculated  from  the  oxygen 
mole  fraction  profile  for  1 ,  2  and  3  bar  cathode  air,  respec¬ 
tively,  as  the  current  density  varies  from  0.05  to  1.2  A  cm-2. 
The  equivalent  oxygen  diffusivities  in  the  cathode  gas  mix¬ 
ture  vary  so  slightly  with  the  current  density  that  they  can 
be  regarded  as  constant.  The  equivalent  oxygen  diffusivities 
were  estimated  as  0.263,  0.117  and  0.075,  respectively,  for 
1 ,  2  and  3  bar  cathode  air. 


3.  Two-dimensional  gas  diffusion  layer  model 

The  two-dimensional  model  for  the  oxygen  transport 
through  the  GDF  takes  the  current  collector  ribs  into  account, 
as  shown  by  the  enlarged  portion  of  the  single  cell  appearing 
in  Fig.  1.  L\  and  L2  denote  the  midline  of  the  GDF  portion 
in  contact  with  the  rib  and  the  channel,  respectively.  A  GDF 
element  bounded  by  L\ ,  L2  has  been  selected  for  analysis  in 
this  two-dimensional  model. 

3.1.  Governing  equation 

The  oxygen  concentration  along  the  catalyst  layer  is  deter¬ 
mined  by  that  in  the  gas  flow  channel  and  by  the  mass  trans¬ 
fer  process  within  the  GDF.  The  two-dimensional  oxygen 
mass  transport  within  the  GDF  can  be  described  by  rewriting 
Eq.  (1 1)  in  vector  form  as: 

No  =  -(e3/2£>tmix)VC  (12) 

where  N0  is  the  oxygen  molar  flux  (a  vector  quantity)  in  the 
GDF.  By  taking  divergence  of  both  sides  of  Eq.  (12)  and 
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by  applying  the  species  conservation  for  oxygen,  we  obtain 
the  Laplace’s  equation  that  governs  the  oxygen  concentration 
distribution  within  the  GDL: 


2  92C  92C 

V  C_  9^+  9?“° 


(13) 


3.2.  Boundary  conditions 


The  boundary  conditions  for  the  GDL  element  under  study 
are  summarized  as: 

3  C  Wr  +  Wc 
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The  catalyst  layer  can  be  regarded  as  an  infinitely  thin 
film  [1,13-15]  located  on  the  left  boundary  of  the  GDL. 
The  rate  of  the  electrochemical  reaction  within  the  catalyst 
layer  can  be  described  using  this  thin  film  model,  yielding  a 
Butler- Volmer  rate  expression.  It  is  then  simplified  to  give  a 
Tafel  type  equation  in  terms  of  the  oxygen  concentration  at 
the  catalyst  layer  as: 
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(14)  3.3.  Solution  procedure 


In  Eq.  (14),  i  is  the  local  current  density,  Av  the  specific 
area  of  the  active  surface,  io  the  reference  exchange  current 
density,  8C  the  thickness  of  the  catalyst  layer,  C(y,  8 d)  the 
oxygen  concentration  at  the  catalyst  layer,  Cref  the  reference 
oxygen  concentration  associated  with  io,  ac  the  cathode  trans¬ 
fer  coefficient,  and  r]c  the  cathode  overpotential. 

In  this  two-dimensional  model,  both  the  local  current  den¬ 
sity  and  the  oxygen  concentration  at  the  catalyst  layer  vary 
withy.  Furthermore,  under  the  steady-state  condition  and  the 
zero  reactant  crossover  assumption,  the  current  is  determined 
by  the  oxygen  diffusion  rate  at  the  GDL/catalyst  layer  inter¬ 
face.  This  indicates  that: 


i  o/9  p  3C 

4^  =  No\z=Sd  =  —£  D0_m\x  — 


(15) 


Combining  Eq.  (14)  with  Eq.  (15),  we  obtain  the  following 
boundary  condition  for  the  GDL/catalyst  layer  interface: 


3  C  _  Avi08c  C(y,8d)  /  acFrjc\ 

&  z=Sd  ~  “4F£3/2 D^.mix  Cref  exp  \~RT~ ) 

=  -KxC(y,Sd )  (16) 


where  the  parameter  K  is  defined  as: 

_  Avi0Sc  &xp(acFr]c/ RT) 

-  4F£3/2DtmlxCref  (  ' 

The  K  value  can  be  evaluated  for  the  given  cathode  over¬ 
potential  r]c  and  the  physical  parameter  and  property  values: 
Av,  io,  8C,  ac,  T,  s,  De0_mix  and  Cref. 

It  is  appropriate  to  assume  symmetrical  boundary  condi¬ 
tions  on  both  the  upper  and  lower  boundaries  of  the  GDL 
element;  hence  the  condition  3C/3y  =  0  is  imposed  on  both 
L\  and  L 2.  The  molar  flow  rate  across  the  GDL/rib  interface 
is  zero,  so  the  3C/3z  =  0  condition  applies  to  that  part  of  the 
boundary.  For  simplicity,  we  set  C  =  Co  on  the  boundary  fac¬ 
ing  the  channel  and  assumed  that  the  cathode  overpotential 
r\ c  is  constant  along  the  catalyst  layer. 


The  two-dimensional  Laplace’s  equation  (Eq.  (13)) 
associated  with  the  corresponding  boundary  conditions 
(Eqs.  (18)-(22))  has  been  discretized  using  the  finite- 
difference  approach  and  solved  using  an  alternating-direction 
explicit  (ADE)  method  [16].  From  the  solution,  we  obtain  the 
oxygen  concentration  distribution  in  the  GDL,  i.e.  C(y,z)-  The 
oxygen  molar  flow  rate  at  the  GDL/catalyst  layer  interface, 
N0\z=$d,  and  the  local  current  density,  i,  are  then  evaluated 
using  Eq.  (15). 

To  indicate  the  performance  of  the  partly-covered  GDL  in 
transporting  oxygen  from  the  gas  flow  channel  to  the  catalyst 
layer,  an  assessment  may  be  made  by  evaluating  the  effec¬ 
tiveness  of  the  GDL,  §g.  It  is  defined  as  the  ratio  of  the  mass 
transport  rate  through  the  GDL  element  (partly  covered  by 
the  rib)  to  the  mass  transport  rate  obtained  without  the  rib. 
Because  the  mass  transport  rate  of  oxygen  is  proportional  to 
the  current  density,  §g  can  be  expressed  in  terms  of  the  current 
densities  as: 

t,  =  '-f  (23) 

where  I  denotes  the  cell  current  density  without  ribs,  and  /ave 
represents  the  average  current  density  with  ribs.  This  was 
evaluated  using  the  following  equation: 


1 

(wr  +  wc)/2 


/ 


(wr+wc)/2 

0 


i(y)dy 


(24) 


4.  Results  and  discussion 

4.1.  Model  validation 

To  validate  the  two-dimensional  GDL  model  associated 
with  the  thin-film  assumption  for  the  catalyst  layer,  the  fuel 
cell  voltage  is  calculated  and  a  comparison  is  made  between 
the  simulation  results  and  experimental  data. 
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- + - Model  with  constant  resistance 

A  Experimental  results  [17] 

Model  with  adjusted  resistance,  Eq.  (26) 


Current  Density  (A/cm2) 

Fig.  5.  Model  fit  for  the  two-dimensional  model  to  the  experimental  data. 


Neglecting  the  overpotential  on  the  anode  side,  the  fuel 
cell  voltage  is  calculated  as: 

^cell  =  Voc  —  0c  —  ve^  (25) 

where  Voc  is  the  open-circuit  voltage  and  9i  is  the  ohmic  re¬ 
sistance  through  the  fuel  cell.  Fig.  5  shows  how  our  model 
fits  with  the  experimental  data,  given  by  Springer  et  al.  [17], 
for  a  PEM  fuel  cell  with  0.3  mm  GDL,  5  [xm  catalyst  layer, 
and  fed  with  H2/air  at  1  bar.  The  physical  parameters  and 
properties  for  the  fit  are  listed  in  Table  2.  A  substantial  dis¬ 
crepancy  between  the  model  fit  with  constant  91  =  0. 14  Q  cm2 
and  the  experimental  data  has  been  found  at  higher  current 
densities.  For  the  limiting  current  density,  however,  a  good 
match  can  be  seen  between  the  model  fit  and  the  experi¬ 
mental  data.  The  discrepancy  at  higher  current  densities  can 
be  attributed  to  the  increased  ohmic  resistance  of  the  dry- 


Table  2 

Physical  parameter  and  property  values  used  in  two-dimensional  model 


Physical  parameter/property 

Value 

Cell  temperature,  T  (K) 

353.15 

Cathode  transfer  coefficient,  ac 

2.0 

Thickness  of  catalyst  layer,  8C  (|jim) 

5.0 

Thickness  of  gas  diffusion  layer,  <$d  (cm) 

0.03 

Width  of  current  collector  rib,  wY  (cm) 

0.15 

Width  of  gas  flow  channel,  wc  (cm) 

0.15 

Reference  exchange  current  density  times  area, 

5  x  10-4 

Avi0  (A  cm  3) 

Reference  oxygen  concentration,  Cref 

4.62  x  1CT6 

(mol  cm  3) 

Oxygen  concentration  at  GDL/channel 

4.62  x  1CT6 

interface,  Cq  (mol  cm  3) 

Porosity  of  gas  diffusion  layer,  s 

0.25 

Cathode  pressure,  Pt  (bar) 

1.0 

Equivalent  oxygen  diffusivity,  D^_mix  (cm2  s_l ) 

0.263 

Nitrogen-oxygen  mole  ratio  in  gas  flow  channel 

3.76 

Ohmic  resistance,  (Q  cm2) 

0.14 

ing  membrane  when  the  fuel  cell  is  operating  at  high  current 
densities,  or  to  the  ohmic  resistance  arising  in  the  catalyst 
layer  that  has  been  neglected  in  the  thin-film  model.  Account¬ 
ing  for  these  effects,  we  have  adjusted  the  ohmic  resistance 
value  with  the  increase  in  current  density  using  the  following 
relation: 


91  =  0.14  + 0.1  (/ave)2  (26) 

As  can  be  seen  from  Fig.  5,  a  good  agreement  has  been 
achieved  for  the  operating  range  up  to  the  limiting  current 
density  of  1 . 1 5  A  cm-2 . 

4.2.  Concentration  distribution  and  current  density 
variation 

The  oxygen  concentration  distribution  within  the  GDL 
element  accompanied  by  the  variation  in  local  current  density 
i  along  the  catalyst  layer  is  shown  in  Fig.  6  for  three  average 
current  densities.  The  same  parameters  listed  in  Table  2  were 
used  for  the  simulations. 

For  the  average  current  density  of  0. 1  A  cm-2,  minor  vari¬ 
ations  in  the  oxygen  concentration  (ranging  from  0.78  Co 
to  Co)  and  local  current  density  i  along  the  catalyst  layer 
(ranging  from  0.087  to  0.107  A  cm-2)  have  been  observed. 
Mathematical  reasoning  has  suggested  that  the  local  current 
density  i  at  y  =  0  be  identical  to  the  current  density  for  the 
PEM  fuel  cell  without  current  collector  ribs  and  hence  §g  has 
been  found  to  be  as  high  as  93%. 

Fig.  6(b)  gives  the  results  for  the  average  current  density  of 
0.5  A  cm-2.  Because  of  the  increase  in  average  current  den¬ 
sity,  a  higher  concentration  gradient  is  required  to  transport 
more  oxygen  to  the  catalyst  layer.  A  substantial  variation  in 
oxygen  concentration  has  been  found  (ranging  from  0.2  Co 
to  Co).  The  local  current  density  for  the  portion  covered  by 
the  rib  is  substantially  lower  than  that  facing  the  channel.  The 
local  current  density  ranges  from  0.22  to  0.7 1  A  cm-2  and  §g 
drops  to  70%. 

As  the  average  current  density  increases  further  to 
1.0  A  cm-2,  from  Fig.  6(c)  the  oxygen  concentration  can 
be  seen  to  have  varied  drastically  and  only  a  trace  of  oxy¬ 
gen  can  be  found  in  the  inner  portion  covered  by  the  rib. 
At  the  catalyst  layer  for  y  >  wc/ 2,  the  local  current  density 
i  decreases  from  1.2  to  0.11  A  cm-2,  which  is  only  about 
1/15  of  the  current  density  at  y  =  0,  and  §g  is  as  low  as 
60%.  This  phenomenon  can  be  attributed  to  the  GDL  being 
so  thin  (8d  =  0.3  mm),  as  compared  to  the  rib  width  (wY  = 
1.5  mm),  that  the  small  cross-sectional  area  of  the  path  re¬ 
stricts  the  mass  transfer  from  the  gas  flow  channel  to  the  y  > 
wc/ 2  portion,  especially  when  the  total  mass  transfer  rate  is 
high. 

The  results  shown  in  Fig.  6  are  based  on  the  assumption 
that  the  cathode  overpotential  r]c  is  constant  along  the  catalyst 
layer.  But  when  the  ohmic  voltage  drop  across  the  membrane 
is  taken  into  account,  the  non-uniformity  in  the  current  den¬ 
sity  might  be  mitigated. 
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4.3.  Influence  of  wr  and  wc  on  GDL  effectiveness 

To  investigate  the  rib  and  channel  size  effects  on  the  oxy¬ 
gen  mass  transport  through  the  GDL,  a  series  of  model  simu¬ 
lations  has  been  conducted  to  evaluate  the  GDL  effectiveness 
with  the  rib  and  channel  widths  being  varied  while  all  of  the 
other  parameters  remained  the  same  as  listed  in  Table  2. 

Fig.  7  shows  the  variation  in  §g  versus  the  cell  current  den¬ 
sity  for  nine  sets  of  wY  and  wc  combinations.  The  circles  at 
the  end  of  each  curve  indicate  that  the  limiting  current  condi¬ 
tions  were  reached  for  the  given  wY  and  wc  combination.  No 
more  data  points  are  expected  to  appear  further  right.  It  is  ob¬ 
vious  that  §g  decreases  with  the  rib  width  and  increases  with 
the  channel  width.  All  of  the  §g  curves  fall  with  the  current 


density,  indicating  that  the  GDL  is  utilized  less  effectively  at 
higher  current  densities.  Because  the  GDL  portion  covered 
by  the  rib  can  still  be  utilized  to  a  certain  extent  even  at  high 
current  densities,  the  §g  value  at  the  limiting  current  density  is 
slightly  higher  than  the  channel  width  fraction  wc/(wY  +  wc). 

4.4.  GDL  thickness  effects  on  cell  performance 

For  the  rib-covered  GDL  portion  (y  >  wc/ 2),  faster  mass 
transport  can  be  achieved  by  increasing  the  GDL  thickness 
because  it  will  enlarge  the  path  for  oxygen  transport  in  the 
y-direction.  However,  as  for  the  channel-facing  portion  (y  < 
wc/ 2),  the  oxygen  mass  transport  rate  decreases  with  the 
increase  in  GDL  thickness.  Thus,  there  may  exist  an  optimal 


(a)  Iave  =  0.1  A/cm2 


(b)  Lve  =  0.5  A/cm2 


Fig.  6.  Distribution  of  oxygen  concentration  within  the  GDL  element  and  the  variation  in  local  current  density  along  the  catalyst  layer  (the  numerical  labels  on 
the  iso-concentration  curves  designate  their  C/Co  values). 
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Current  Density  i  (A/cm2)  Dimensionless  Position  ( zlSd ) 

(c)  lave  =  1.0  A/cm2 
Fig.  6.  ( Continued ). 


thickness  for  the  GDL  under  a  given  set  of  wv  and  wc ,  and 
this  possibility  will  be  further  explored  below. 

Fig.  8  shows  the  calculated  performance  curves  for  a  single 
PEM  fuel  cell  of  various  GDL  thicknesses  with  all  of  the  other 
parameters  being  identical  to  those  listed  in  Table  2.  As  can  be 
seen,  the  overall  performance  falls  with  the  increase  in  GDL 
thickness.  A  very  slight  discrepancy  in  cell  voltage  occurs 
at  low  current  densities,  but  substantial  deviation  appears  at 
mid  and  high  current  densities.  The  limiting  current  density 
decreases  as  the  thickness  increases. 

The  current  density  at  Vceii  =  0.6  V  (a  typical  operating 
voltage  of  PEM  fuel  cells)  has  been  singled  out  for  quantita¬ 


tive  investigations.  The  current  density  versus  GDL  thickness 
is  plotted  in  the  inset  in  Fig.  8  together  with  the  other  two  sets 
of  wv  and  wc.  The  relationship  among  the  three  curves  in  the 
inset  reflects  their  relative  ranking  in  terms  of  §g,  as  depicted 
in  Fig.  7.  The  curves  in  the  inset  fall  with  the  increase  of 
GDL  thickness  and  no  optimal  thickness  for  the  GDL  has 
been  found. 

We  should  bear  in  mind  that,  in  addition  to  thickness, 
the  porosity  of  the  GDL  also  influences  the  effective  cross- 
sectional  area  of  the  oxygen  transport  path  in  the  y-direction. 
Fig.  9  shows  the  cell’s  current  density  (again  at  Vcen  =  0.6  V) 
versus  GDL  thickness  of  various  GDL  porosities,  with  both 


Fig.  7.  GDL  effectiveness  vs.  current  density  for  various  rib  and  channel  widths. 
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Fig.  8.  Calculated  performance  curves  for  a  single  PEM  fuel  cell  at  various  GDL  thicknesses. 


the  rib  and  the  channel  measured  1.5  mm  in  width.  It  can 
be  seen  from  the  curves  that  the  performance  of  a  PEM  fuel 
cell  has  substantially  enhanced  with  an  increase  in  the  GDL 
porosity.  As  the  porosity  increases  to  0.5  or  above,  a  max¬ 
imum  current  density  emerges  and  the  maxima  are  desig¬ 
nated  on  each  curve  by  circles.  The  curves  in  Fig.  9  indicate 
that,  from  the  mass  transfer  point  of  view,  the  GDL  thickness 
should  be  as  thin  as  possible  if  the  GDL  porosity  is  low.  How¬ 
ever,  an  optimal  GDL  thickness  exists  when  a  high-porosity 
GDL  is  used. 

For  electrical  conductivity,  material  strength  and  ease  of 
fabrication  considerations,  the  GDL  generally  has  a  porosity 
substantially  lower  than  1.0.  Under  such  a  design  restriction, 
an  increase  of  the  oxygen  mass  transport  in  the  y-direction 
may  offer  the  possibility  of  enhancing  the  overall  effective- 


Fig.  9.  Cell  current  density  (at  0.6  V)  vs.  GDL  thickness  for  various  GDL 
porosities. 


ness  of  the  GDL.  This  may  be  carried  out  by  using  a  GDL 
that  has  a  directional  preference  in  transporting  oxygen.  For 
instance,  yarns  of  carbon  cloth  parallel  to  the  channel  flow 
direction  are  substantially  thinner  than  that  of  perpendicular 
yarns,  and  a  sparser  arrangement  can  be  made  for  yams  to  be 
perpendicular  to  the  channel  flow. 


5.  Conclusions 

We  have  developed  a  two-dimensional  model  that  simpli¬ 
fies  the  descriptions  of  kinetic  and  mass  transfer  processes  in 
PEM  fuel  cell  cathode.  This  model  is  validated  with  experi¬ 
mental  data  from  the  literature  and  is  used  to  investigate  the 
mass  transfer  problem  within  the  GDL  under  the  influence  of 
current  collector  ribs. 

The  GDL  can  be  used  effectively  at  low  current  densi¬ 
ties;  however,  its  effectiveness  falls  with  increasing  current 
density.  At  high  current  densities,  the  effectiveness  of  the 
GDL  is  only  slightly  higher  than  the  channel  width  fraction 
wc/(wT  +  wc).  The  mass  transfer  is  slow  in  the  GDL  portion 
covered  by  the  ribs  and  only  a  trace  of  oxygen  can  be  found 
in  the  inner  part  when  the  current  density  is  high. 

From  the  mass  transfer  point  of  view,  the  GDL  thickness 
should  be  as  thin  as  possible  if  the  GDL  porosity  is  low. 
However,  an  optimal  thickness  for  the  GDL  exists  when  a 
high-porosity  GDL  is  used.  The  existence  of  such  an  optimal 
GDL  thickness  implies  that  a  GDL  with  directional  prefer¬ 
ence  in  mass  transfer  can  be  used  to  improve  PEM  fuel  cell 
performance. 
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